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Abstract 

This paper proposes an approach based on Bacterial 
Foraging Algorithm (BFA) for optimal placement of 
Distribution Generations (DGs). The optimal sitting and 
sizing of distributed generation is formulated as a multi- 
objective function including the network power losses 
minimization and voltage profile improvement. BFA 
algorithm has been employed to optimize the problem. To 
demonstrate the effectiveness of the proposed approach, it 
has been applied to 33 and 69 bus systems. The results of the 
proposed approach are compared with those of Genetic 
Algorithm (GA) and also Shuffled Frog Leaping Algorithm 
(SFLA) to reveal its strong performance. 
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Introduction 

Distributed generation is expected to play an 
increasing role in emerging electrical power systems. 
Studies have predicted that DG will be a meaningful 
percentage of all new generations going on lines. It is 
predicted that they are about 20% of the new 
generations being installed [1]. Optimization 
techniques should be employed for deregulation of 
the power industry, by applying the best allocation 
of the distributed generations (DGs). The 
advancement in technology and the demand of the 
customers for cheap and reliable electric power has 
led to an increasing interest in distributed 
generation. The issues related to reliability and 
maintenance has impeded the penetration of DG 
resources in distribution systems [2-6]. The DG units 
might bring different benefits such as: network 
investment deferral [7], active loss reduction [8-10], 
voltage profile improvement [22], and reliability 
improvement [12-14]. The benefits of DG units highly 



depend on the size and location of them in the network. 
Many methods have been proposed in the literature to 
find the optimal location and size of DG units in the 
network which have considered various technical 
aspects such as: voltage limits, feeder capacity limits 
and number of installed DG units [15]. In [16], a GA 
method is proposed to minimize the distribution 
systems active power flow. In [17], a combination of 
GA and simulated annealing meta-heuristic methods is 
used to solve optimal DG power output. A TS method 
to size the DG optimally, as well as the reactive sources 
within the distribution system, is presented in [18]. A 
multi objective function is proposed to determine the 
optimal locations to place DGs in distribution system 
to minimize power loss of the system and enhance 
reliability improvement and voltage profile [19]. [20] 
Considered an analytical expression to calculate the 
optimal size and an effective methodology to identify 
the corresponding optimum location for DG placement 
for minimizing the total power losses in primary 
distribution systems. The selection of the best places 
for installation and the preferable size of the DG units 
in large distribution systems is a complex 
combinatorial optimization problem. In [21], 
Lagrangian based approaches are used to determine 
optimal locations for placing DG, considering 
economic limits and stability limits. In [22-24], the 
authors succeeded in merging both the DG location 
and size in one optimization problem. The main factors 
included in the optimization problem were investment 
cost, operation cost, network configuration, active and 
reactive power costs, heat and power requirements, 
and voltage profile and system losses. Several methods 
have been adopted to solve such an optimization 
problem. Some of them rely on conventional 
optimization methods and others use artificial 
intelligent-based optimization methods. In some 
researches, the optimum location and size of a single 
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DG unit is determined [22], whereas in others the 
optimum locations and sizes of multiple DG units are 
determined [23]. In [25], a mixed integer linear 
program was formulated to solve the optimization 
problem. The objective was to optimally determine the 
DG plant mix on a network section. However, that 
required dealing with the power system 
approximately as a linear system which is not the real 
case. In [26], a tabu search (TS)-based method was 
proposed to find the optimal solution of their problem, 
but the TS is known to be time-consuming algorithm in 
addition to its ability to be trapped in a local minimum. 
In [22], a particle swarm optimization (PSO) algorithm 
was introduced to determine the optimum size and 
location of a single DG unit to minimize the real power 
losses of the system. The problem was formulated as 
one of the constrained mixed integer non-linear 
programming with the location being discrete and the 
size being continuous. However, the real power loss of 
the system was the only aspect considered in this work 
while trying to optimally find the size of only one DG 
unit to be placed. In [23], different scenarios were 
suggested for optimum distribution planning. One of 
these scenarios was to place multiple DG units at 
certain locations pre-determined by the electric utility 
distribution companies (DISCOs) aiming to improve 
their profiles and minimize the investment risk. In this 
paper, the problem of optimal DG location and sizing 
in distribution systems is formulated as a multi- 
objective optimization problem and BFA is used to 
solve this problem. The optimal DG location and sizing 
problem is converted to an optimization problem with 
the multi-objective function including the network 
power losses and better voltage regulation. The 
effectiveness of the proposed BFA is tested on 33 and 
69 bus systems in comparison with the GA and SFLA 
through some performance indices. Results show that 
the proposed method provides the correct answers 
with high accuracy in the initial iterations in 
comparison with other methods and is superior to 
them. 

Allocation ofDG 

The distribution systems are usually regulated through 
tap changing at substation transformers and by the use 
of voltage regulators and capacitors on the feeders. 
This form of voltage regulation assumes power flows 
circulating from the substation to the loads. DG 
introduces meshed power flows that may interfere 
with the traditional used regulation practices. Since the 
control of voltage regulation is usually based on radial 



power flows, the inappropriate DG allocation can 
cause low or over-voltages in the network. On the 
other hand, the installation of DG can have positive 
impacts on the distribution system by enabling 
reactive compensation for a voltage control, reducing 
the losses, contributing for frequency regulation and 
acting as spinning reserve in main system fault cases. 

Heuristic Optimization Methods 

Genetic Algorithm (GA) 

It is well known that GAs work are according to 
the mechanism of natural selection. Stronger 
individuals are likely to be the winners in a 
competitive environment. In practical applications, 
each individual is codified into a chromosome 
consisting of genes, that representing a characteristic of 
one individual. For identification parameters of a 
model, parameters are regarded as the genes of a 
chromosome, and a positive value, generally is known 
as the fitness value. It is used to reflect the degree of 
goodness of the chromosome. Typically, a 
chromosome is structured by a string of values in 
binary form, which the mutation operator can 
operate on any one of the bits, and the crossover 
operator can operate on any boundary of each two 
bits in the string [27]. Since, the parameters in our 
problem are real numbers, a real coded GA is used, 
in which the chromosome is defined as an array of real 
numbers with the mutation and crossover operators. 
Here, the mutation can change the value of a real 
number randomly, and the crossover can take place 
only at the boundary of two real numbers. 

Shuffled Frog Leaping Algorithm (SFLA) 

SFLA is a meta heuristic optimization method that 
mimics the memetic evolution of a group of frogs 
when seeking for the location having maximum 
amount of available food. The algorithm contains 
elements of local search and global information 
exchange [28]. The SFLA involves a population of 
possible solutions defined by a set of virtual frogs 
partitioned into subsets referred to as memeplexes. 
Within each memeplex, the individual frog holds ideas 
that can be influenced by the ideas of other frogs, and 
the ideas can evolve through a process of memetic 
evolution. Using a particle swarm optimization like 
method the SLFA simultaneously performs an 
independent local search in each memeplex. To 
ensure of global exploration, after a defined number 
of memeplex evolution steps (i.e. local search 
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iterations), the virtual frogs are shuffled and 
reorganized into new memeplexes in a technique 
similar to that used in the shuffled complex 
evolution algorithm. In addition, to provide the 
opportunity for random generation of improved 
information, if the local search cannot find better 
solutions as random, virtual frogs are generated and 
substituted in the population. The local searches and 
the shuffling processes continue until defined 
convergence criteria are satisfied. The SFLA is 
described in details as follows. First, an initial 
population of N f rogs P = { X l , X 2 , . . . , X N } is cr eated 
as random.For S-dimensional problems (S variables), 
the position of a frog i' h in the search space is 
represented as X i = [x x , x 2 ,...,x is ] T . 

Afterwards, the frogs are sorted in a descending order 
according to their fitness. Then, the entire population 
is divided into m memeplexes, that each containing n 
frogs (i.e. N=mxn), in such a way that the first frog 
goes to the first memeplex, the second frog goes to the 
second memeplex, the m th frog goes to the m ,h 
memeplex, and the (m + l)' h frog goes back to the first 
memeplex, etc. 

Let M k is the set of frogs in the k th memeplex, this 

dividing process can be described by the following 
expression: 

M k = {X k+m{l _ X) ^P\l<k<n},(l<k<m). (1) 

Within each memeplex, the frogs with the best and the 
worst fitness are identified as X h and X w , 

respectively. Also, the frog with the global best fitness 
is identified as X g . During memeplex evolution, the 
worst frog X w leaps toward the best frog X h . 

According to the original frog leaping rule, the 
position of the worst frog is updated as follows: 

D=r.(X b -X w ), (2) 
X w (new) = X w + D,(\\D\\<D m J, (3) 

Where r is a random number between and 1; and 
D max is the maximum allowed change of frog's 

position in one jump. Fig.4 demonstrates the original 
frog leaping rule. If this leaping produces a better 
solution, it replaces the worst frog. Otherwise, the 
calculations in (2) and (3) are repeated with regard to 
the global best frog. If in this case, no improvement 
becomes possible, the worst frog is deleted and a new 



randomly generated frog is replaced. The calculations 
continue for a predefined number of memetic 
evolutionary steps within each memeplex, and then 
the whole population is mixed together in the shuffling 
process. The local evolution and global shuffling is 
continued until a convergence criterion is satisfied. 
Usually, the convergence criteria can be defined as 
follows: 

i. The relative change in the fitness of the best frog 
within a number of consecutive shuffling iterations is 
less than a pre-specified tolerance; 

ii. The maximum user-specified number shuffling 
iterations is reached. 

The SFLA will stop when one of the above criteria is 
achieved first. 

Bacterial Foraging Algorithm (BFA) 

Natural selection tends to eliminate animals with poor 
foraging strategies and favor the propagation of genes 
of those animals that have successful foraging 
strategies since they are more likely to enjoy 
reproductive success. After many generations, poor 
foraging strategies are either eliminated or shaped into 
good ones. The E. coli bacteria that are present in our 
intestines also undergo a foraging strategy. The control 
system of these bacteria that dictates how foraging 
should proceed can be subdivided into four sections 
namely Chemotaxis, Swarming, Reproduction and 
Elimination and Dispersal [29]. 

1) Chemotaxis 

This process is achieved through swimming and 
tumbling via Flagella. Depending upon the rotation of 
Flagella in each bacterium, it decides whether it should 
move in a predefined direction (swimming) or 
altogether in different directions (tumbling), in the 
entire lifetime. To represent a tumble, a unit length 
random direction, say cf>, is generated; this will be used 
to define the direction of movement after a tumble in a 
tumble, the position of the ith bacterium is updated as: 

& (j + lkJ) = i (j,k,l) + CxZ(/> (4) 

Where 6\j,k,l) is the position of the ith bacterium at 
the jth chemotactic step of the kth reproduction loop in 
the 1* elimination-dispersal event, C is the size of the 
step taken in the random direction specified by the 
tumble, /-(/> is the angle of the direction which is 
randomly generated in the range of [0, 2ri) . 
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2) Cell-to-cell Communications 

E-coli bacterium has a specific sensing, actuation and 
decision-making mechanism. As each bacterium 
moves, it releases attractant to signal other bacteria to 
swarm towards it. Meanwhile, each bacterium releases 
repellent to warn other bacteria to keep a safe distance 
from it. BFA simulates this social behavior by 
representing the combined cell-to-cell attraction and 
repelling effect as: 

J cc (9 i ( j,k,l),9( j,k,l)= {.J'cc(9 l ,99) 
t=\ 



-- I 
t=\ 

s 

- I 
t=t 



~ d attract ex P( '^attract M0 m~9 m) ) 
m=l 



~d re pellant ex P(~ w repeUant Z(9'm ~9 t m ) 2 ) 
m=\ 



(5) 



Where J cc (0',0) is the cost function value to be 

added to the actual cost function to be minimized to 
present a time varying cost function. 'S' is the total 
number of bacteria and 'P' the number of parameters to 
be optimized which are present in each bacterium. 

dattract, COattract ,drepellant ,GJrepellant are different coefficients 

that are to be chosen properly. 
3) Reproduction 

In BFA, a fixed total number of reproduction steps, Nre 
is given. Only the first half of populations survive in 
each reproduce step a surviving bacterium splits into 
two identical ones, which occupy the same position in 
the environment as the one in previous step. Thus the 
population of bacteria keeps constant in each 
chemotactic step. After Nc chemotactic steps, the fitness 
values for the ith bacterium in the chemotactic loop are 
accumulated and calculated by: 



. Nc+l • 

j health = Z j J ( j,k,l) 
7=1 



(6) 



Where j 'health presents the health of the ith bacterium, 

the smaller the J health IS, the healthier the bacterium is. 
To simulate the reproduction character in nature and 
to accelerate the swarming speed, all the bacteria are 
sorted according to their health values in an ascending 
order and each of the first Sr ( Sr = S / 2 ), for 
convenience S is assumed to be a positive even integer) 
bacteria splits into two bacteria. The characters 
including location and step length of the mother 
bacterium are reproduced to the children bacteria. 
Through this selection process the remaining Sr 
unhealthier bacteria are eliminated and discarded. To 



simplify the algorithm, the number of the bacteria 
keeps constant in the whole process. 

4) Elimination-dispersal 

For the purpose of improving the global search ability, 
elimination-dispersal event is defined after Nre steps of 
reproduction. The bacteria is eliminated and dispersed 
to random positions in the optimization domain 
according to the probability Ped. This elimination- 
dispersal event helps the bacterium avoid being 
trapped into local optima. The number of the event is 
denoted as Ned. The flowchart of the BFA has been 
depicted in Fig. 1. 

Problem Formulation 

The optimal placement and sizing problems of 
distributed generation are formulated as a multi- 
objective function including the network power losses, 
better voltage regulation and improve the voltage 
stability. The goal is to converge these three objective 
functions into one, using the penalty coefficient. 

Fitness Function 

The objective function is calculated as follows [24]: 

OF = Min ((A + kj 2 + fe 2 / 3 ) + ft Z ieWDC [max(V ni - V™ x ,0) + 
maxVnimax- Vni, 0+fi2i£NmaxSni-Snimax, (7) 

Where fi, ti and f3 are the objective functions to the real 
power losses, improve voltage profile and for 
improving voltage stability index, respectively which 
are expressed in components as follow [30]: 

fl = * ( p gm ~ Pdni ~ V mi V ni Y m cos( 8 mi - 8 ni + 9 ni )) (8) 

i=2 



fi— T.(y n i -y rated y 

ni=l 



1 



(S/(«0) 



ni=2,3, . . .iin 



(9) 
(10) 



Wheres 



SKn 2 ) =\V mi I 4 -4[PnAnW ni +Q ni {ni)X m ]\V mi I 2 -4[,P„, (11) 

It is very important to identify weak buses for nodes 
with minimum voltage stability index that are prone 
to voltage instability. Investigating the voltage 
stability index behavior demonstrate that the buses 
which experiencing large voltage drops are weak and 
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within the context of remedial actions. So, it makes 
sense to act on controls that will improve the voltage 
magnitudes at weak buses. 

1) Power Flow 

The power flow equations to be satisfied for each bus 
as follow [24]: 

N 

Pgni ~ P dni ~ Vni £ V nj Y nj cos( S ni - 8 nj - 9 nj ) = (12) 
7=1 

N 

Qgni ~ Qdni ~ ^ni ^-^nj^nj s ^ n ( &ni ~ &nj ~ @nj ) = ® 

1) Voltage Profile 

The voltage should be kept within standard limits at 
each bus as follow [24]: 



tt nun y tt max 
v ni ^ v ni ^ v ni 



(14) 



2) DG Units 



It is necessary for the DG units to be constrained 
capacity between the maximum and the minimum 
levels as follow [24]: 



p mm p p n 



3) Thermal constraint 



(15) 



Thermal limit of distribution lines for the network 
should not be exceeded [24]: 



(16) 



■5*^" Stop ^) 





FIGURE 1 FLOWCHART OF BFA 

Case Study 

The proposed BFA scheme for optimal placement and 
sizing of DG has been tested on several power systems. 
Here, the test results for two different distribution 
systems are presented to demonstrate its effectiveness 
and the results are investigated. The values of the 
objective functions before the installation of DGs have 
been shown in Table 2. 

TABLE 2 OBIECTIVE FUNCTION VALUE OF THE SYSTEMS 
BEFORE DG INSTALLATION 



System 



Objective function value 



Fl 



F2 



F3 



33-Bus 
69-Bus 



0.2109 
0.2217 



0.3141 
0.2197 



1.4907 
1.4509 



1) Coding and Decoding of the Solution 

The coded variables are, therefore, the candidate locals 
for DG installation and the respective when applying 
genetic algorithms to optimize the DG allocation and 
sizing problem, an important aspect is the coding of 
the potential solutions. In a general way, a potential 
solution is a configuration of the network with the DG 
units size of the units. In that sense, each chromosome 
is coded through a vector whose dimension is the 
number of candidate points and whose content is the 
size of the installed units. The chromosome coding 
used in the proposed methodology is of the integer 
type {0, 1, 2, 3,}. In the DG allocation process, number 
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The adjustable value of BFA, GA, SFLA algorithms for simulation are given in table 1. 

TABLE 1 ADJUSTED PARAMETERS OF BFA, GA, SFLA ALGORITHMS 



BFA parameters 




GA parameters 


number of bacteria 


6 


population size 50 


dimension of optimization 
process 


6 


number of optimization 
variables 


number of swimming step 


8 


Noise band 0.5*[3 3 3 1 1 1] 


number of chemotaxis step 


30 


mutation rate 0.8 


number of reproduction 


5 


cross over rate 0.2 


number of elimination and 




Creation of initial 

, . rand(50,6) 
population 


dispersal 


5 


probability index of 
elimination and dispersal 


0.35 




normalized distribution 


Rand (6 , 6) 




step size for each 
bacterium 


0.05 




elimination-dispersal loop 
index 







reprodution loop index 







chemotaxis loop index 








means no DG unit, while the other numbers represent 
sequentially a DG unit from a list of candidates 
provided for each candidate installation place. The 
decoding of the GA solution is based on the same idea, 
translating a number in the chromosome into the 
respective DG unit. 

2) 33 Bus Radial Distribution Systems 

At first the proposed technique is applied on a radial 
system using the total load 3.72 MW, 2.3 MVar, 33 
bus and 32 branches as it is shown in Fig2. The real 
power and reactive power losses in the system are 
210.998 (kW) and 143 (kVar) respectively. The 
voltage profile before installation DG and after 
optimally placing DGs, and voltage stability index 
are given in Fig. 3 and 4, respectively. The minimum 
fitness value evaluating process has been shown in 
Fig. 5. Table 3 demonstrates the results for optimal 
placement and sizing problems of distributed 
generation. 



3) 69 Bus Radial Distribution Systems 

The other system employed to evaluate the proposed 
method is the 69 bus radial distribution system that has 
the total load of 3.80 MW and 2.69 MVar and it is 
demonstrated in Fig. 6. For more details, the voltage 
profile and voltage stability index are illustrated in 
Fig. 9 and Fig. 8. The minimum fitness value evaluating 
process DG location and capacities are represented in 
Fig 11 and Table 4. 
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FIGURE 3 VOLTAGE PROFILE OF 33-BUS 
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FIGURE 4 VOLTAGE STABILITY INDEX OF 33-BUS 
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FIGURE 6 THE 69 BUS RADIAL DISTRIBUTION SYSTEM 



The results presented in table 2 and table 3, shows that 
the BFA is effective for optimal placement and sizing of 
DG. In addition, the running time for the BFA was 
faster in comparison with the others. As it obvious 
from the Fig. 5 and fig. 7, the voltage profiles at all 
nodes for radial transmission systems have improved. 
Fig. 5 and fig. 8 depict Convergence profile of each bus 
in 33 and 69 bus distribution systems. It can be seen 
that the convergence of BFA algorithm is faster 
than SFLA and GA algorithm and this is because 
BFA algorithm provides the correct answers with 
high accuracy in the initial iterations which makes 
the responding time of this algorithm extremely fast. 
The voltage stability for two systems has been shown 
in Fig 4 and 8. As it obvious from the figures, the 
voltage stability after installing DGs has been 
improved. 

4) Stability and Voltage Regulation 

Voltage stability index in bus 18 from the first system 
and bus 61 from the second were low before DG 
installation. This could cause instability in the 
networks in the presence of disturbances. After DG 
installation, the three methods showed major 
improvements, Figs. 4 and 8. Besides, the Figs. 3 and 7 
are demonstrating that the voltage regulation index is 
at highest for BFA and for the GA is the lowest. 

Conclusion 

In this paper, the BFA is presented for optimal 
placement and sizing of DG units. The proposed BFA 
algorithm for optimal placement of DG is easy in 
performance without additional computational 
complexity. The capability of the proposed approach is 
tested on 33 and 69 bus systems to minimize the losses, 
increase the voltage stability and improve the voltage 
profile. The simulation results show that the BFA 
yields has better convergence characteristics and 
performance compared with the other heuristic 
methods such as GA and SFLA. Also the proposed 
technique exhibits a higher capability in finding 
optimum solutions by taking into account the active 
power and reactive power losses for objective function. 
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TABLE 3 PERFORMANCE ANALYSIS OF THE 33-BUS SYSTEM AFTER DC INSTALLATION 



Objective function value 



Method 

XII TIT 133 

rl rl rj 


Bus.no 


DG 
size(MW) 


BFA 0.0921 0.0022 1.0247 


25 


1.1770 




13 


0.7771 




31 


±-±yyj 


SFLA 0.0940 0.0028 1.0236 


31 


1.1944 




12 


0.8894 




5 


1.0664 


GA 0.1005 0.0058 1.0302 


28 


0.6632 




15 


0.8177 




30 


1.1297 


TABLE 4 PERFORMANCE ANALYSIS OF THE 69-BUS SYSTEM AFTER DG INSTALLATION 


Objective function value 


Method 

Fl F2 F3 


Bus.no 


DG 
size(MW) 


BFA 0.0921 0.0022 1.0247 


25 


1.1770 




13 


0.7771 




31 


1.1993 


SFLA 0.0940 0.0028 1.0236 


31 


1.1944 




12 


0.8894 




5 


1.0664 


GA 0.1005 0.0058 1.0302 


28 


0.6632 




15 


0.8177 




30 


1.1297 




FIGURE 7 VOLTAGE PROFILE OF 69-BUS FIGURE 8 VOLTAGE STABILITY INDEX OF 69-BUS 
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FIGURE 9 VARIATIONS OF OBJECTIVE FUNCTION 

Appendix 

rin total number of buses in the given RDS 

ni receiving bus number (ni =2,3, ... ,n) 

mi bus number that sending power to bus m (rm 
= ru =1) 

i branch number that fed bus m 

Ndg total number of DG 

Pgni active power output of the generator at bus m 

Qgni reactive power output of the generator at bus m 

Pdni active power demand at bus m 

Qdni reactive power demand at bus m 

Pni(ni) total real power load fed through bus m 

Qni(ni) total reactive power load through bus m 

P mm Minimum active power of DG at bus ni 

5 ni 1 
D max , t->.^t 

r ni Maximum active power of DG at bus ni 

Vni voltage of bus m 

Vmi voltage of bus mi 

V ni mm Minimum voltage at bus m 

V ni Maximum voltage at bus m 

I S m max I Maximum apparent power at bus m 
Yni admittance between bus ni and bus mi 

9 ni Phase angle of Yi= Y ni Z ni 

S ni Phase angle of voltage at bus m V ni = V ni Z. 5 ni ) 



8 mj Phase angle of voltage at bus mi 

Vreted rated voltage (1 p.u.) 

SI(m) voltage stability index of node m. 

(ni =2,3,... ,n) 

pi penalty coefficient, 0.32 

p2 penalty coefficient, 0.3 

n penalty coefficient (ki = 0.6) 

r2 penalty coefficient (k2= 0.35) 

fi network real power losses (pu) 

f2 network voltage profile (pu) 

f3 network voltage stability index (pu) 
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